Accepted for publication in ApJ, Draft version September 9, 2011 

Preprint typeset using L£Tf<]X style cmulatcapj v. 11/10/09 



OBSERVATIONS OF ENERGETIC HIGH MAGNETIC FIELD PULSARS 
WITH THE FERMI LARGE AREA TELESCOPE 

D. Parent 1 , M. Kerr 2 , P. R. den Hartog 2 , M. G. Baring 3 , M. E. DeCesar 4,5 , C. M. Espinoza , E. V. Gotthelf 7 , 
A. K. Harding 4 , S. Johnston 11 , V. M. Kaspi 14 , M. Livingstone 14 , R. W. Romani 2 , B. W. Stappers 6 , K. Watters 2 , 
P. Weltevrede b , A. A. Abdo 1 , M. Burgay 13 , F. Camilo 7 , H. A. Craig 2 , P. C. C. Freire 10 , F. Giordano 8,9 , 
L. Guillemot 10 , G. Hobbs 11 , M. Keith 11 , M. Kramer 6 ' 10 , A. G. Lyne 6 , R. N. Manchester 11 , A. Noutsos 10 , 

A. Possenti 13 , D. A. Smith 12 

Accepted for publication in ApJ, Draft version September 9, 2011 

ABSTRACT 

We report the detection of 7-ray pulsations from the high-magnetic-field rotation-powered pul- 
sar PSR J1119— 6127 using data from the Fermi Large Area Telescope. The 7-ray light curve of 
PSR J1119— 6127 shows a single, wide peak offset from the radio peak by 0.43 ± 0.02 in phase. Spec- 
tral analysis suggests a power law of index 1.0 ± 0.3+Q2 with an energy cut-off at 0.8 ± 0.2+05 GeV. 
The first uncertainty is statistical and the second is systematic. We discuss the emission models of 
PSR J1119— 6127 and demonstrate that despite the object's high surface magnetic field — near that 
of magnetars — the field strength and structure in the 7-ray emitting zone are apparently similar to 
those of typical young pulsars. Additionally, we present upper limits on the 7-ray pulsed emission 
for the magnetically active PSR J1846— 0258 in the supernova remnant Kesteven 75 and two other 
energetic high-5 pulsars, PSRs J1718— 3718 and J1734— 3333. We explore possible explanations for 
the non-detection of these three objects, including peculiarities in their emission geometry. 
Subject headings: gamma rays: stars — pulsars: general — pulsars: individual (PSR J1119— 6127, 
PSR J1718-3718, PSR J1734-3333, PSR J1846-0258) 



1. INTRODUCTION 

Magnetic fields in neutron stars are thought to origi- 
nate mainly from the fossil fields of the progenitor stars 
or from fi eld generatio n during supernova core collapse 
(see, e.g.. iSpruitl [2008). These strong fields can man- 
ifest themselves as 'classical' rotation-powered pulsars 
and magnetars, also known as anomalous X-ray pulsars 
(AXPs) and soft 7-ray repeaters (SGRs). While pul- 
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sars emit steady, beamed electromagnetic radiation from 
radio up to high energies ultimately via the loss of ro- 
tational kinetic energy (E = 4ir 2 IP/P 3 ergs -1 ), mag- 
netars dissipate their extremely high surface magnetic 
fields (B s ~ 10 14 — 10 15 G) in luminous X-ray emission — 
variable on multiple time scales — generally exceeding the 
power derived from losing rotational kineti c energy (see 
I Woods fc Thompsod200l[Mereg hetti 20Ql for reviews) . 

The bulk of pulsar magnetic fields inferred from spin 
parameters are clustered around 10 12 G. Nevertheless, 
a few have estimated fields near or surpassing the 
quantum-critical field B CI = rn 2 c 3 /eh = 4.4 x 10 13 G, 
at which the electron cyclotron energy is equal to its 
rest mass. This small class of strong-field pulsars pro- 
vides an opportunity to constrain emission mechanisms 
at high energy and to explore the interesting in t erface 
between pulsars and magnetars. ICamilo et all pOOOh 
discovered PSR J1814— 1744, the first pulsar breaking 
the B CI barrier with B s = 5.5 x 10 13 G and E = 4.7 x 
10 32 ergs -1 , and PSR J1119-6127 with a slightly weaker 
field. The latter object recently exhibited an unusual 
glitch recovery during which its radio emission shifted 
from its typical profile, showing an "intermittent" peak 
and "Rotating- Radio- Tr ansient-like" c haracteristics (see 
lMcLaughlin"et~aT1 [20061 for RRATs). iWeltevrede etaLI 
(|2010al) demonstrated that this type of behavior could 
be related to reconfigurations of the magnetic field in the 
magnetosphere. Moreover, PSR J1846-0258, the pulsar 
in the supernova remnant (SNR) Kesteven 75, exhibited 
the first known magnetic activity in a rotation-powered 
pulsar with a scries of 5 SGR-likc X-ray bursts and an 
X-ray brightening ( both pulsed and DC emission) lasting 
about two months (jGavriil et al J 120081 ) . A large spin- up 
glitch with an unusual recovery coinciding with the on- 
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set of its magnetar-like behavior has also been rep orted 
(jKuiper fc Hermsenl [200l ILivingstone et al.1 I2010D . In 
this context, the connection between high- 5 pulsars and 
magnctars is strengthened by the discovery of a low- 
B SGR which shows ma gnetar-like activ ity despite its 
canonical magnetic field (jRea et al.ll2010f ). 

In the high energy field for high-5 pulsars, the lower 
energy spectral cut-off between 10 and 30 MeV observed 
for PSR B1509-58 in the Energetic Gamma-Ray Ex- 
periment Telescope (EGRET) er a was attributed to 
the higher surfac e magnetic field ([Harding et al.l 119971 : 
iThompsonl I2008T ). However, the recent discovery of 
GeV pulsations from the pulsar in the supernova rem- 
nant CTA1 with a comparable magnetic field and no 
radio counterpart sh owed that the em ission mechanism 
is more complicated (jAbdo et al.ll2008h . These observa- 
tions clearly encourage closer examination of these ob- 
jects to help assess if this class has unique properties 
from a possible overlap with magnetars. 

In this paper we present 7-ray observations, acquired 
by the Large Area Telescope (LAT) aboard the Fermi 
Gamma-ray Space Telescope, of energetic high-i? pulsars 
(B > 4 x 10 13 G) for which the spin-down energy loss rate 
is above 10 33 ergs _1 . We report the detection of pulsed 
7-ray emission from PSR J1119— 6127 and discuss this 
detection and possible emission mechanism scenarios. 
We also present an upper limit on the pulsed emissions of 
the three other candidates including PSR J1718— 3718. 
PSR J1734-3333, and PSR J1846-0258. Table Q] lists 
measured and derived parameters for the selected pul- 
sars, ordered by their surface field strengths B s . 

PSR J1119— 6127 was discovered in the Parkes multi- 
beam pulsar survey with a period of 408 ms and a 
large period derivative of 4.0 x 10~ 12 s s _1 ([Camilo et al.l 
120001) . Subsequently, it was detected by XMM- Newton 
in the 0.5 - 2.0 kcV range sho wing a single, narrow 
pulse aligned with the radio peak ([Gonzalez et al.1120051 ). 
Besides its large surface dipole magnetic field B s of 
4.1 x 10 13 G, the measurement of timing parameters in- 
fers a very young pulsar (t c = P/2P) of 1.6 kyr with 
a high spin-down power of E = 2.3 x 10 36 ergs -1 . 
PSR J1119— 6127 is associated with the ^15' diame- 
ter shell-type SN R G292.2-0.5 observed in both radio 
and X -ray bands ([Crawford et al.ll2001t IPivovaroff et al.l 
l2001t ). as well as a compact and faint X-ray pulsar 
wind n ebula (PWN) observed by Chandra X-ray Obser- 
vator y ([Gonzalez fc Safi-Harbll2003l : ISafi-Harb fc Kumarl 
120081 ) which is spati ally coincident with t he Te V source 
HESS J1119-614 (|Diannati-Atai et al.1 [20091 ). The 
SNR/PSR J1119-6127 system is located in the Galactic 
plane 8.4 ± 0.4 kpc away, based on neutral hydrogen ab - 
sorption of X-rays from the SNR (jCaswell et all 120041 ) . 
However, using th e NE2001 model of th e Galactic elec- 
tron distribution ([Cordes fc Laziol 120021 ) . we derived a 
distance of 16. 7^^ kpc, which places the pulsar beyond 
the Sagittarius arm. This discrepancy is probably due 
to the line-of-sight [l = 292.151, b ='-0.537), which is 
nearly aligned with Sagittarius arm, making distance de- 
termination more uncertain than usual. In this work, we 
use the distance d = 8.4 kpc. 

PSR J1718-3718 is a young pulsar with a large spin 
period of 3.3 s and a period derivative of P = 1.6 x 
10 _12 ss _1 , implying a spin-down luminosity of 1.6 x 



10 ergs , the lowest of the four selected pulsars. How- 
ever, it has the second highest surface magnetic field 
strength (Bs = 7.4 x 10 13 G) among the known ra- 
dio pulsars. It was discovered in the Parkes Multi- 
beam Survey ([Hobbs et al.l l2004| ) not far from the di- 
rection of the Galactic Center (I = 349.85, b = 0.22), 
and observed as a point-like source by Chandra X-ray 
Observatory with a thermal spectrum resembling that 
of the transient AXP XTE J1810-197 in quiescence 
([Kaspi fc McLaughlin! 120051 ). Recently, deeper Chandra 
X-ray Observatory observations have shown pulsations 
in the soft X-ray band ( 0.8 - 2keV) and better con- 
strain the source spectrum ([Zhu et al.ll2011l) . Its distance 
based on a dispersion measure (DM) of 373pccm~ 3 is 
d = 4.5 ±0.5 kpc. 

PSR J1734-3333 is a radio pulsar fo und in the Parkes 
Multib eam Survey ([Morris et all 12002ft with B s = 5.2 x 
10 13 G. lOlausen et "all ([20101 ) report the probable X-ray 
detection of the pulsar using XMM-Newton observations. 
The NE2001 model assigns a distance of 6.1^ 1 Qkpc 
which, with the projected position (I = 354.82, b = 
—0.43), implies the pulsar is located in the Near 3 kpc 
Arm toward the Galactic Center. 

PSR J1846-0258 was disco vered in the Rossi X-ray 
Timing Explorer data f RXTE: iGotthelf et al.|[2000l) with 
B s = 4.9 x 10 13 G. Its spin period of 324ms and its spin- 
down rate of P = 7.1 x 10 _12 ss _1 indicate this pul- 
sar is one of the youngest known to date (t c = 723 yr). 
No radio pulse d emission has been d etected yet despite 
deep searches ([Archibald et al.1 [20081 . The distance of 
the associated SNR Kesteven 75 and the pulsar is still 
in debate. Published estimates have ranged from 6.6 kpc 
to 21 kpc (Caswell et al. 1975, Milne 1 979, Becker fc 
Helfan d 1984). A more recent study by iLeahy fc Tianl 
(|2008f ) estimates the distance between 5.1 kpc and 7.5 kpc 
based on HI a nd 13 CO velocity measurements, while 
ISu et all (p009f ) in an detailed study of molecular-cloud 
velocities from 12 CO measurements determine the dis- 
tance as 10.6^°'gkpc. Based on these references, we 
adopt in this work the distance of 7.9 ± 2.8 kpc. 

2. OBSERVATIONS 

2.1. Timing Observations 

Among the pulsars considered here, three have spin- 
down powers above 10 34 ergs -1 (see Table [1} and there- 
fore have been followed by several observatories in the 
context of the pulsar timing campaign for Fermi, which 
monitors the mos t energetic radio and X-ray pulsars 
(|Smith et al.ll2008f ). This extensive program, begun in 
early 2007 and lasting through the Fermi LAT mission, 
provides accurate measurements of the pulsar rotation 
parameters used to assign phases to 7-ray photons. 

PSRs J1119-6127 and J1734-3333 are observed ap- 
proximately monthly with the 64-m Parkes radio tele- 
scope in Australia, at a frequency of 1.4 GHz (and oc- 
casional observations at 0.7 and 3.1GHz). The Jodrell 
Bank Observatory also monitors PSR J1734— 3333 on 
a weekly basis with the 76-m Lovell telescope, using 
a 64 MHz band centered at 1404 MHz connected to an 
analog filterbank. For these pulsars a total of 94 and 
78 pulse times of arrival (TOAs) were recorded from 
mid-2008 to 2010, overlapping the F ermi LAT observa - 
tions. The TEMP02 timing package (jHobbs et al.ll2006ft 
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was used to build the timing solutions from the TOAs, 
which have been fitted to the pulsar rotation frequencies 
and their derivatives. For PSR Jl 119-6127, the fit in- 
cludes 16 harmonically related sinusoids, using the "FIT- 
WAVES" option in the TEMP02 package, to flatten the 
timing noise. The post-fit rms is 2.1 and 30.7ms, or 
0.5% and 2.6% of the pulsar phase, for PSRs J1119-6127 
and J1734— 3333 respectively. Full det ails of the observ- 
ing and dat a analysis can be fo und in iWeltevrede et all 
(|2010bt) and lHobbs et~aT1 (pOOl . 

PSR J1718-3718, due to its relatively low spin- 
down power, is not a target for the Fermi campaign. 
We used a timing solution constructed as outlined in 
iManchester fe Hob bs (2010). Timing observations were 
made between 2009 January 24 and 2009 June 19 us- 
ing the Parkcs radio telescope with the 10-cm receiver, 
which has a 1024-MHz bandwidth centered at 3100 MHz, 
and the Parkcs digital filtcrbank system. The resulting 
timing solution based on 23 TOAs has a post-fit rms of 
3.4 ms. 

Finally, X-ray timing observations for PSR J1846-0258 
were made weekly using the Proportional Counter Ar- 
ray (PCA; Uahoda et aLlll996t ) on RXTE, and were col- 
lected in "Good Xenon" mode. Photons arriving be- 
tween 2008 June 5 and 2010 April 6 were extracted 
from the first xenon detection layer in the energy range 
2 - 20kcV 1 and folded with the previously published 
ephemeris (jLivingstone et aI]|2010D . The resulting pulse 
profiles were cross-correlated with a high significance 
template, producing a single TOA for each observa- 
tion. The TOAs were fitted to a timing model using 
the pulsar timing software package T EMPCp^l. Further 
detail s of t he analysis can be f ound in iLivingstone et aT] 
(|2006h and ILivingstone et al.1 ()2010l ). Note that the de- 
tailed description of the timing analysis and results for 
PSR J1119— 6127 using the PCA instrument are given in 

The timing parameters used in this work will be made 
available from the Fermi Science Support CcntciP^I. 

2.2. Gamma-ray Observations 

The LAT aboard Fermi is an electron-positron pair 
conversion telescope and went into orbit on 2008 June 11 
(jAtwood et all l2009t lAbdo et al.1 12009Q . The telescope 
covers the 20 MeV to > 300 GeV energy range with 
good sensitivity and localization performance (an effec- 
tive area ~ 8000 cm 2 on-axis above 1 GeV and an angular 
resolution 9qs ~ 0.6° at 1 GeV for events in the front sec- 
tion of the tracker). The LAT timing is derived from a 
GPS clock on the spacecraft, and 7-rays are hardware 
time- stamped to an a ccuracy significantly better than 
1 fjs (|Abdo et al.H2009f) . 

3. ANALYSIS AND RESULTS 
3.1. Data Selection 

For each pulsar, we selected LAT data collected be- 
tween 2008 August 4 (MJD 54682) when Fermi began 
scanning-mode operation^3 and the end of the ephemeris 

° http:/ /www. atnf.csiro.au/research/pulsar/tempo 

16 http: / /fermi. gsfc.nasa.gov / ssc / data/access/lat /ephems 

17 Except for J1718-3718 for which the timing solution starts at 
54855 MJD. 



validity range (17-29 months). We used "diffuse" class 
events (highest probability of being 7-ray photons) un- 
der the P6_V3 instrument response function (IRFs), and 
excluded events with zenith angles > 100° to reject at- 
mospheric 7-rays from the Earth's limb. In addition, we 
also have excluded for the spectral analysis time intervals 
when the region of interest (ROI) intersects the Earth's 
limb. The events were analyzed usin g t he standard 
software package Science Tools-09-21-00E\ (ST) for the 
Fermi LAT data analysis and photon phases were calcu- 
lated using the "ferrm- plugin" available in the TEMP02 
pulsar timing software (|Rav et al.ll2011| ). 

Multiple scattering at low photon energy dominates 
the LAT's angular resolution. To approximate as much 
as possible the instrument Point Spread Function (PSF), 
we selected events for the light curves with an energy- 
dependent angular radius cut centered on the stud- 
ied pulsar. This selection is given by < Qq${E) > = 
[(5.12°) 2 x (100MeV/E) ls + (0.07 ) 2 ] 1 / 2 which approx- 
imates a 68% containment angle according to the IRFs. 

We searched for 7-ray pulsations using the bin- 
independent H-test (|de Jager fc B iisching 201(|), and se- 
lecting photons passing the energy-dependent cut for dif- 
ferent maximum angular radii (from 0.5 to 5 degrees) 
around the pulsar position and energy bands. This trun- 
cates the point-spread function at low energies and de- 
creases the number of background events. After taking 
the number of trials into account, we detected a pulsed 
7-ray signal from PSR J1119— 6127 with a significance 
above 5 a (our detection threshold). For the other pul- 
sars we found significances below 3 tr as discussed in § 13.51 

3.2. Gamma-ray Study of PSR Jl 119- 6127 
3.2.1. Light Curves 

Figured] presents the 20 bin 7-ray phase histograms in 
three energy ranges, along with the radio profile and the 
X-ray profile observed by XMM-Newton (bottom pan- 
els). The top panel shows the profile with optimized 
signal-to-noise ratio, corresponding to an event selection 
above 0.5 GeV within 0?5 of the timing position and pass- 
ing the energy-dependent cut. For this energy band, 
the H-test gives a value of 107 corresponding to a pul- 
sation significance of 9 a. The dashed line shows the 
background level (76 counts per bin) estimated using an 
annular ring centered on the radio position, during the 
off-pulse window (</> < 0.10 and <p > 0.65), and with inner 
and outer radii of 0?5 and 1?5, respectively. No signif- 
icant pulsed signal was detected below 0.5 GeV, despite 
searches using different apertures. The light curve con- 
sists of one single, wide peak observed between 0.10 and 
0.65 in phase, which defines our on-pulse window. The 
light curves between 0.5 — 1 GeV and above 1 GeV within 
the same aperture have similar pulse profiles. We fitted 
the unbinned 7-ray data above 0.5 GeV with a Gaussian 
function. The peak is offset in phase from the radio peak 
by 5 = 0.43 ± 0.02 with a FWHM of 0.18 ± 0.03. The 
bias due to the DM uncertainty in extrapolating the ra- 
dio TOA to infinite frequency is negligible. 

3.2.2. Spectrum 

18 http: / /fermi. gsfc.nasa.gov/ssc/data/analysis/scitools/overview. html 
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Fig. 1. — Phase-aligned radio, X-ray, and 7-ray profiles of 
PSR J1119-6127. Bottom panel: 1.4 GHz radio profile observed 
with the Parkes telescope. Second panel from bottom: X-ray pro- 
file in the 0.5-2.0 ke V ranges observed by XMM (adapted from 
Gonzalez ct al. 2005). Top panels: 7-ray profiles obtained with 
the Fermi Large Area Telescope in different energy bands. Two 
rotations are shown. The dashed line shows the background level 
estimated from a surrounding ring. 



A phase-averaged spectr um was obtained with a 
binned likelihood analysis (jMattox et al.l 119961 ) of the 
LAT data selected from 20° x 20° region centered on 
the pulsar position, using the Fermi Science Tool "gt- 
like". The Galactic diffuse emission was modeled using 
the glLiem_v02 map cube, while the extragalactic emis- 
sion and residual instrument backgrounds were modeled 
jointly by the isotropic component isotropic_iem_v02. 
These two models and an explanation of their develop- 
ment are available from the FSSC. In addition, all point 
sources within 15° of the pulsar found in an updated 
version of the LAT 1FGL source catalog (jAbdo et al.l 
I2010bl) . using 18 months of data, were included in the 
model. Sources were modeled with a power law spec- 
trum, except for pulsars, for which a power law with an 
expon ential cut-off was used (see EqQ] and lAbdo et"aL1 
l2010dl) . Sources more than 5° from the pulsar were as- 
signed fixed spectra taken from the source catalog. Spec- 
tral parameters for sources within 5° of the pulsar were 
left free for the analysis. We fitted the spectrum of 
PSR J1119-6127 above 100 MeV using an exponentially 
cut-off power-law of the form 



dN 
~dE 



= N E 



-r. 



exp 



E 
E~ c 



cm -2 s -1 MeV" 1 , 



(1) 



where Nq is the differential flux (in units of 
phcrn -2 s _1 MeV -1 ), T the photon index, and E c the 
cut-off energy. The parameter /3, which determines the 
steepness of the exponential cut-off, was fixed to 1. The 
Fermi LAT pulsars are generally well-described by a 
simple exponential model, /3 = 1, a shap e predicted 
by ou ter magnetosphere emission models (jAbdo et al.l 
I2010dn . The energy at which the normalization factor 
No is defined is 1 GeV. The best-fit values are listed in 
Table Q] as well as both the photon flux Fioo and energy 
flux Gioo above 100 MeV, where the first errors are sta- 
tistical and the second are systematic. The systematic 
uncertainties were estimated by applying the same fit- 
ting procedures described above and comparing results 
using bracketing IRFs where the effective area has been 
perturbed by ± 10% at 0.1 GeV, ± 5% near 0.5 GeV, 
and ± 20% at 10 GeV with linear interpolations in the 
logarithm of the energy for intermediate energies. We 
also modeled the pulsar with a simple power-law spec- 
trum (/3 = 0), and by leaving the f3 parameter free in 
Equation [TJ In the first case, the exponentially cut-off 
power law model is preferr ed at the 6 a level according 
to the likelihood ratio test (|Mattox et al.|[l996l ). while in 
the second case the extra free parameter did not improve 
the quality of the fit. The results were all cross-checked 
using an alternate analysis tool developed by the LAT 
team. Figure [2] shows both the phase-averaged spectral 
fit between 0.1 and 10 GeV (solid lines) with /3 = 1, and 
the spectral points derived from likelihood fits to each in- 
dividual energy band in which it was assumed the pulsar 
had a power-law spectrum. 

To search for unpulscd emission from a compact nebula 
or magnetospheric emission from the pulsar, we searched 
in the off-pulse region (tp < 0.10 and <fi > 0.65) for a 
point source in the energy band 0.1 — 100 GeV at the 
radio pulsar position. Using a power-law spectrum we 
found a la signal with an energy flux above 0.1 GeV of 
(2.0±0.3)xl0 -11 ergs -1 which represents ~ 30% of the 
phase-averaged emission. There are three possibilities 
to explain this signal: emission from a PWN or SNR, 
magnetospheric pulsar emission, or unmodclcd structure 
in the Galactic diffuse emission. To determine its ori- 
gin, we fitted the signal in three energy bands (0.1 - 
1 GeV, 1 - 10 GeV, 10 - 100 GeV) and in the whole en- 
ergy range (0.1 - 100 GeV) using a pulsar shape spectrum 
(see eq. [T]) . The signal is only significant below 1 GeV 
and the exponentially cut-off power law model is pre- 
ferred over a simple power-law model at the 5 a level, 
with a spectral index of 0.2 ± 1.0 and an energy cut- 
off of 0.3 ± 0.2 GeV. An exponentially cut-off spectrum, 
whic h has been observed in the off-pulse of a few pul- 
sars (|Ackermann et aLll2011l ) , suggests a magnetospheric 
origin. However, the model for the Galactic diffuse emis- 
sion is imprecise in this complex region, and such fea- 
tures may also be associated with unmodeled structure 
in the Galactic diffuse emission. In particular, the ap- 
parent low energy cut-off may be a result of the large 
PSF at low energies (6» 68 ~ 3° at 100 MeV for events 
in the front section of the tracker). If the emission is 
in fact an inadequately-modeled diffuse component, the 
phase-averaged 7-ray flux for PSR J1119— 6127 would 
be reduced by 30%. However, because the emission is 
only present in a narrow band at low energy, we expect 
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the phase-averaged spectral shape to remain largely un- 
changed. 

Finally, the pulsar associated to the source 
1FGL J1119.4— 6127c does not show any variabil- 
ity in the 1FGL catalog and an internal 18-month 
LAT source list (|Abdo et all I2010b( ). But in view of 
the recent evidence for changes in the magnetospheric 
configuration of PSR J1119— 6127, a careful search for 
changes in the spectrum and pulse shape is warranted if 
a glitch occurs. 




Energy(GeV) 

Fig. 2. — Phase-averaged spectral energy distribution of 
PSR J1119— 6127 in 7-rays obtained with the Fermi Large Area 
Telescope. Plotted points arc from likelihood fits to individual en- 
ergy bands with > 4 cr detection above background for two degrees 
of freedom, otherwise an upper limit arrow is shown. The solid 
black line shows the maximum likelihood fit to a power law with 
exponential cut-off (Eq. [TJ. The dashed lines are ± 1 a uncertain- 
ties on the fit parameters. 



3.3. X-ray Study of PSR Jl 119- 6127 

We also searched for pulsed X-ray emission from 
PSR J1119— 6127 to compare to our Fermi results. The 
pulsar was monitored regularly as part of an RXTE mon- 
itoring campaign^ to look for magn etar-like behavior a s 
discovered from PSR J1846-0258 (jGavriil et al.1 120081 ). 
Data were collected using the PCA, an array of five colli- 
matcd xenon/methane multi-anode proportional counter 
units (PCUs) sensitive to incoming photons in the 2 - 
60 keV energy range, with a total effective area of approx- 
imately 6500 cm 2 and a field of view of - 1° FWHM. We 
used all archival XTE data available in HEASARC pub- 
lic archivf£3 on the pulsar that spanned the valid interval 
time of the radio ephemeris. All data were taken in high 
time resolution Good- Xenon mode (>100 fis). These 
data sets were processed using the standard reduction 
and analysis tools. After extracting time intervals con- 
taminated with PCA breakdown events, the data were 
merged into a single barycenter corrected FITS photon 
file. We selected counts from the first layer and energy 
range up to 20 keV and folded the data with the same 
ephemeris used to analyze the Fermi LAT data (see §3.1). 
No significant signal was found in the XTE data range 
(MJD 54890 - 55170). This is not surprising since the X - 
ray pulsed emission detected by iGonzalez et al.l ()2005l ) 
shows strong energy dependence and appears to be con- 
fined to the soft X-ray band below 2 keV. We then folded 

19 PI: F. P. Gavriil 

20 http:/ /heasarc. gsfc.nasa.gov/docs/archive. html 



data from energy bands restricted to the lowest energy 
PI channels to test all cumulative channels up to 5 keV. 
However, again no signal was detected. 

3.4. Emission Models of PSR Jl 119- 6127 

The first year of pulsar studies with Fermi has shown 
that the 7-ray emission mechanism for you ng pulsars is 
likely to be outer magnetospheric in origin ([Abdo et al.l 
I2010dfl . Such emission is predicted by s everal differ- 
ent m od els. The Oute r Gap models (OG, iCheng et al.1 
[l986llbl : 1 omanil IT996I ) place the emitting region on a 
set of open zone field lines paralleling the closed zone 
boundary; particles accelerate and radiate between the 
'null charge' surface and the light cylinder, so that only 
one magnetic pole is visible in each hemisphere and the 
hollow cone produces the typical double 7-ray pulse. An 
alternative picture is the two-pole caustic model (TPC, 
IDvks fc Rudakl 12001 which uses a similar set of field 
lines, but posits acceleration starting near the star sur- 
face; both poles are visible in both hemispheres and a 
double pulse is formed by truncating the emission zone 
at ~ 0.8 r^c so that only the trailing edges of the two 
cones are present in the light curv e. More physical pic- 
tures , e.g. the slot-gap model (SG. lMuslimov fc Harding! 
I2001 . may have radiation from a wide range of altitudes 
and so may be approximated by either of these two ge- 
ometries. 

We can compare the predictions o f the OG and TPC 
pictur es following the procedure in iRomani fc Wattersl 
(|2010( ). For these calculations we assume a co-rotating 
vacuum-like field ('pseudo-force-free'), model the field 

for a gap width w = (lO^ergs^ 1 /^) 172 01 0.02 for 
PSR J1119— 6127, and compute the light curve for each 
magnetic inclination angle a and viewing angle C- I n 
Figure [3] we compare the match between the model light 
curve and the dat a in the (a, () plane usi ng the \3 
weighting defined in IRomani fc Wattersl (|2010l ): dark col- 
ors arc better fits. For PSR J1119— 6127, as for many 
young Fermi pulsars, the radio emission has high lin- 
ear polarization and has been fit to the rotating vec- 
tor model ( RVM) to constrain the g eometrical angles a 
and ( (e.g... iCrawford fc Keiml l200ll ) . Even better con- 
straints can be obtained from the u nusual double-peaked 
mode found by iWeltevrede et all (|2010aD in post-glitch 
timing observations. We fit the RVM model to these 
data and correct for the R VM sign convention problem 
([Everett fc Weis berg 20*011) to obtain the allowed range 
(green contours) in Figure [3J As is often the case, these 
fits primarily constrain /3 = £ — a; the 0.5, 1.5 and 2.5 er 
contours from these fits are shown. 

The idealized RVM model assumes a point dipole and 
does not account for the radio emission altitude or the 
magnctosphcre closed zone. If the radio emission arises 
from finite altitude, the phase and shape of the pulse 
maximum and the radio polarization sweep are shifted 
(jBlaskiewicz et al.l Il99lt IDvks fc Harding] l200l . We 
have used our numerical models of the magnetospheric 
field structure to fit the polarization sweep for finite al- 
titude radio emission, requiring that the observed radio 
flux arise from open field lines. The best fits arc found 
for altitudes ~ 0.1 r^c! we show the contours of the al- 
lowed regions in the left (OG) panel of Figure [3] While 
the allowed a, ( range is quite similar to that of the 
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Fig. 3. — Pulsar geometry and emission modeling for PSR J1119— 6127 for the 'Outer Gap' (left panel) and the 'Two Pole Causti c ' (right 
panel) models in the magnetic inclination — viewing angle (a, Q plane. Green contours show the RVM fit to the Weltcvrcdc ct al. (2010a) 
radio polarization data. For the left panel we also show the polarization fit for finite altitude, open zone radio emission (r = 0.09 r^c 
magenta contours). Contours are at 1.5x, 2.5x, and 3.5x the minimum value of the reduced x 2 = 0.85. The background color scale gives 
the X3 statistic fit to the observed > 500 MeV pulse profile. The color scales in the panels are the same, with dark colors representing 
better fits. Preferred models lie along the diagonal polarization fit band. 



RVM fit, the phase constraints provide important addi- 
tional information. In Figure [?] we show the observed 
and model light curves for the radio emission for a and £ 
at the best values allowed by the 7-ray pulse shape and 
polarization constraints in Figure [3] The lower panel 
shows the Parkes radio light curve in the double-peaked 
mode. The corresponding positional angle (PA) sweep 
of the polarization vector is shown in the upper panel 
(right scale). The phasing is set by the fit (green dotted 
line) to the maximum rate of the PA sweep, d\& / d0 max 
(top panel). This fit is very good for the inferred angles 
(x 2 /DoF=l.l) and, in turn, determines the phase of the 
closest approach of the surface magnetic dipole axis to 
the Earth line-of-sight. This is denoted by the pulsar 
phase <j>B = 0. Note that the phase of the observed pulse 
centroid (I ma x, lower panel) is offset from the dipole axis. 
The intensity and PA offse ts are in the same s ense a s the 
analytic approximation bvlBlaskicwi cz et al.l (11991 1 ). but 
field line sweep back distortion ( Dvks fc Hardind I2004D 
makes the offsets small er, implying a higher emis sion al- 
titu de. Note also tha t IWeltevrede et alj (|2010a[ ), using 
the IBlaskiewicz et ail (|1991f ) approximations, derived a 
lower emission height. Thus to fit the observed pulse 
width into the open zone they argued for small a. With 
the numerical modeling of the polarization sweep and 
pulse offset we derive larger heights and find that the 
observed sweep and pulse profile are well matched for 
angles consistent with the 7-ray data, while having the 
radio emission arise from the open zone. 

We can now compare with the observed Fermi 7-ray 
light curve (Figure 01 upper panel). Adding the con- 
straints supplied by the polarization data restricts the 
acceptable region. For the OG model the best 7-ray 
fits consistent with the polarization constraints are at 



a = 125 - 130°, C = 140 - 150°. For the TPC picture 
best fits in the allowed range are also near a = 125°, 
£ = 145°, however this model gives a poorer pulse profile 
(%3 = 8.1 v.s 6.0) fit than the OG case (Figure©. The 
Fermi light curve is referenced to the radio pulse peak (in 
the single-peaked mode) and the observed radio-peak/7- 
pulse offset is S = 0.43. The 7-ray pulse models are com- 
puted relative to the magnetic axis (4>b = 0); we find 
that the best-fit OG light curve (solid line), in fact, has 
phase 4>b = 0.39 for altitudes of w 0.1 Tlc, correspond- 
ing to S = 0.43. The model pulse is single, with little 
off-pulse emission. Using the phase and viewing angles 
determined by the radio observations, the best-fit TPC 
model has two peaks and a large unpulsed component 
(Figure 4). For these models the flux correction factors 
are / n = 1.50 (OG) and f n = 0.95 (TPC). 

Note that these fits assume the background level in- 
ferred from the surrounding ring. If, in contrast, the 
true minimum of the magnctosphcric flux is the baseline 
level at <j> = 0.65 — 0.1, then we find that the OG light 
curves fit much better than the TPC case (x3 = 2-0 vs. 
11.3); the best fit a and £ are nearly unchanged. 

We conclude that, in the context of these simple ge- 
ometrical models, the high altitude component conven- 
tionally associated with the OG model is the preferred 
origin of the observed 7-ray pulse, with an excellent 
match to the observed pulse phase and light curve shape. 
However, as for other single-peaked 7-ray pulsars it is 
difficult to exclude other model scenarios. In particular, 
some versions of the slot gap models that allow emission 
at larger altitudes may provide a good fit to the light 
curve; we will need higher sensitivity to search for the 
fainter low altitude off-pulse emission predicted to oc- 
cur in these models. Also, more physical realizations of 
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Pulsar Phase (0 B ) 

Fig. 4.— Light curves for PSR J1119-6127. The bottom panel 
shows the Parkes radio light curve in the two peaked (post-glitch) 
mode, which provides the best model constraints. The correspond- 
ing radio polarization position angle data are shown (right scale) 
in the upper panel. The model Pulsar Phase (</>b) is referenced 
to the closest approach of the magnetic axis to the Earth line-of- 
sight, as fit from the polarization sweep. The sweep rate maximum 
(green arrow) and pulse profile offsets (magenta arrow) are shown, 
with good matches to the observed radio data for an altitude of 
r = 0.09 tlc- The upper panel shows the LAT pulse profile (left 
scale) and model OG (solid line) and TPC (dashed line) profiles. 
These are best-lit profiles (geometric angles in the legend) and the 
phase is referenced to the radio-determined phase of the magnetic 
axis. The 7-ray background, shown by the dotted line, was esti- 
mated using an annular ring centered on the radio position with 
inner and outer radii of 0.5° and 1.5° respectively, during the off- 
pulse region. 



the pulsar magnctosphcre, such as th e numerical force- 
free 'S eparatrix Layer' (SL) model of iBai fc~S pitkovskv 
(|2010f ) may produce acceptable light curve fits. In sum- 
mary, we find that the light curve of Jl 119— 6127 can 
be well fit by conventional 7-ray pulsar models but only 
when emission from relatively close to the light cylinder 
dominates. Despite the object's high surface magnetic 
field, the field strength and structure in the 7-ray emit- 
ting zone are apparently similar to those of more typical 
young pulsars. 

3.5. Gamma-ray Upper Limits for Other High Magnetic 
Field Pulsars 

To calculate an upper limit on pulsed flux, we sim- 
ulated a point source at the position of the pulsar us- 
ing the Science Tool gtobssim and the P6.V3.DIFFUSE 
IRFs. We used the actual pointing history of the space- 
craft to accurately account for the changes in rocking 
profile during the mission. Such profile changes affect 
the exposure a given source receives and can appreciably 
alter sensitivity. The time intervals of the simulations 
are limited to the range of validity of the timing solution 
for each pulsar, an integration time of about 73, 82, and 
87 weeks for the pulsars J1718-3718, J1734-3333, and 



J1846— 0258. Since there is no detection, pulsed or other- 
wise, of these three pulsars, we adopted a spectral shape 
for the simulated source typical of the shape observed 
for other LAT pulsars, namely a power law with expo- 
nential cut-off (see Eq. [J), with F = 1.5, E c = 3GeV, 
and j3 = 1. We note that higher values of E c and lower 
values of r — i.e., sources with a greater fraction of their 
emission at higher energies, where the LAT has superior 
angular resolution and the Galactic diffuse emission is 
less intense — are easier to detect with both pulsed and 
unpulscd methods. In addition to the candidate pulsar, 
we simulated events from the Galactic (glLiem_v02) and 
isotropic (isotropic_iem_v02) diffuse backgrounds. 

Next we modeled the pulsar's light curve as a sin- 
gle Gaussian peak and assigned to each simulated "pul- 
sar" photon a phase drawn at random from this light 
curve. To the diffuse background photons we assigned 
a phase drawn at random from a uniform distribution. 
To calculate the significance of pulsation, we selected 
photons from a 5 x 5 grid of circular apertures con- 
structed such that all included photons had a recon- 
structed position within 6 deg of the point source and 
a reconstructed energy above a threshold Eth, with 
5 linearly spaced between 0?5 and 1?0, and E t h uni- 
formly spaced in logarithmic energy between 100 MeV 
and 1 Ge Y. For each photon set, we c alculated the H-tcst 
statistic (|de Jager fc Bii sching 2010), and the pulsed sig- 
nificance was determined as the chance probability of the 
maximum observed test statistic multiplied by a "trials 
factor" of 25. This combination of pulsation test statistic 
and grid search yields good sensitivity to pulsations for 
a source of unknown spectrum and unknown light curve 
in an arbitrary background. 

To determine the upper limit on pulsed emission, we 
simulated an ensemble of 20 sources at a series of integral 
fluxes, and for each ensemble we determined the distribu- 
tion of pulsation significances via the above method. The 
upper limit was then simply the flux at which 68% of the 
sources had a pulsed significance above a given thresh- 
old, here taken to be 5c, or a chance probability less than 
5.8 x 10 -7 . We note there is an approximately linear re- 
lationship between the flux threshold derived and the 
significance threshold, so the results reported here can 
be scaled to less stringent detection criteria. Finally, we 
perform this exercise for a series of light curves by vary- 
ing the width of the gaussian peak. The results, shown 
in Figure [5j encompass both the sharp peaks typical of 
7-ray pulsar emission and broader sinusoidal peaks. 

Finally, we verified the accuracy of the technique by 
performing a similar exercise with actual LAT data in 
place of the simulated background. The resulting limits 
agree closely with those derived from simulation, indicat- 
ing the impact of neglecting the contributions of other 
point sources to the background is unimportant in deriv- 
ing the pulsed detection thresholds. 

3.6. Gamma-ray Luminosity and Efficiency 

As argued by lHarding fe Muslimovl (|2002[ ). pulsars 
that produce electron-positron pairs through curvature 
radiation have their primary acceleration limited by the 
effect of screening of the electric field. This property 
implies that the high-energy luminosity 
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Fig. 5. — The detection threshold for pulsed emission assuming 
a 100% pulsed fraction from three pulsars as a function of peak 
width. The width is characterized by the FWHM of single gaussian 
peak of the assumed light curve, and the flux reported is that 
for which 68% of the simulated sources have a pulsed significance 
(calculated according to the method described in the main text) 
greater than 5o\ The integral flux reported is for photon energies 
above 100 MeV. 



is proportional to E 1 / 2 , and that 7-ray efficiency 77 7 = 
Ly/E increases with decreasing spin-down power. In 
Eq. [2j d is the pulsar distance, fn is the geometrical 
correction factor depending on the inclination and view- 
ing angles a, £ (see i j3.4p . and G100 is the energy flux 
measured above 0.1 GeV. Figure [6] shows for the studied 
pulsars, as well as the normal 7-ray pulsars reported in 
the first Fermi-LAT Pulsar Catalog (|Abdo et al.ll2010d[ ) 
and iTheureau et al.l (|2011| ). the L 1 as a function of E. 
The reported values are based on the pulsar distances 
discussed in this paper, and under the assumption that 
fn = 1 in order to compare the results with the pulsar 
catalog. 

For PSR J1119— 6127, summing the distance and 
G100 uncertainties in quadrature yields L 1 = (54 ± 
6 ± 8)xl0 34 /si ergs -1 . This places the pulsar be- 
tween the heuristic constant voltage line = 10 33 x 
(E'/IO 33 ) 1 / 2 ergs -1 represented by the dot-dashed line 
and the line for 100% conversion efficiency (L 7 = E). 
This value confirms that in the range 10 35 ergs -1 < E < 



1Q 36.5 ej-gg-i^ 



Lj seems flat. 



For the emission models 
discussed in £|3.4i the derived flux correction factors are 
fa = 1.50 (OG) and jfa = 0.95 (TPC), and hence the 
luminosity cited above may be underestimated for the 
OG model. 

For PSR J1846— 0258, we derived from the upper limit 
on pulsed flux a luminosity of (26 ±19) x 10 34 /n ergs -1 , 
using a distance of 7.9 ± 2.8 kpc and assuming a Gaus- 
sian peak with a FWHM of 0.2. The upper limit overlaps 
the line for small distances (~ 5 kpc) , which is more 
constraining if PSR J1846— 0258 behaves as a rotation- 
powered pulsar. However for large distances in the direc- 
tion of the far side of the Sagittarius arm, the upper limit 
is well above the constant line, matching the luminosity 
of PSR J1420-6048. 



For PSRs J1718-3718 and J1734-3333, the derived 
upper limits on the pulsed luminosity are well above 
100% efficiency, and are thus unconstraining. However, 
these values may be biased by the adopted distances, 
which are based on the DM. 

4. DISCUSSION 

The physics close to the surface of a neutron star with 
a magnetic field comparable to the quantum-critical field 
(B a = 4.4 x 10 13 G) requires a full QED (quantum elec- 
trodynamics) treatment (see lHarding fe La3 l2006L for a 
review) , and we outline a few of the salient features here. 
With this strong field, the cyclotron energy (huj c ) equals 
the rest mass of the electron and is much larger than 
the Coulomb energy. Particles flow primarily along field 
lines, and particle momentum transverse to the field lines 
is quantized in Landau levels. Resonant absorption and 
emission into / from these states is a primary feature of the 
radiation spectra. Most importantly, the strong field en- 
ables otherwise forbidden processes such as single photon 
pair creation 7 — > e + e - and photon splitting 7 — > 77 to 
proceed, and in fact become quite probable. Intense pair 
creation is usually invoke d as a probable seed for signif- 
icant radio emission fe.g. lSturrocklll971l ). Yet these two 
processes can act to suppress the overall yield of created 
electron-positron pairs, in turn possib ly leading to radio 
quiet high- B pulsars and magnctars (jBaring fe Harding! 
119981 ): photon splitting, because it can operate below 
the 7 — > e + e~ threshold of 2m e c 2 , and pair creation, 
because when B > 6 x 10 12 G, it gene rates pairs in the 
ground state (jBaring fe Harding! 120011 ). thereby inhibit- 
ing cascading and subsequent pair creation. The pre- 
cise level of reduction of the pair yield by these physi- 
cal mechanisms is contingent upon a number of factors, 
especially the photon emission/attenuation locale, and 
is therefore somewhat uncer tain, as discussed at length 
bv lBaring fe Harding! (|2001| ). Mor eover, radio emission 
from sever al high-jj pul sars dCamilo et al.l 120001) and 
magnetars (|Camilo et all 120061 l2007t iLevin et alj|2010D 
has been detected, a direct indication that the radio sig- 
nal is not totally suppressed. Perhaps it originates from 
a higher altitude where the magnetic field is weaker and 
pair suppression is limited. 

For 7-ray emission originating near the surface of 
the neutron star, the influence of the high magnetic 
field appears as a super-exponential cut-off in the 7- 
ray spectr um if one-photon pair-creat ion attenuation 
dominates (jDaughertv fe Hardind 119821 ) , while a more 
radual cut-off is exp ected if photon-splitting dominates 
Harding et alJll997j ). These low altitude processes can- 
not modify 7-ray emission originating in the outer mag- 
netospherc. As we argue below, the similarity of the ob- 
servational characteristics of PSR J1119— 6127 to other 
"normal" 7-ray pulsars implies any connection between 
the high surface magnetic field and high-altitude 7-ray 
emission is weak. This is underpinned by the exponen- 
tial character of the turnover in Figure [5] As discussed 
for other Fermi LAT pulsars, the maximum energy e max 
of the pulsations provides a lower limit to the altitude 
of 7-ray emission, since it must lie below the 7-B mag- 
netic absorption threshold. Assuming a dipole field, this 
provides an estimate of the minimum emission altitude 
r > (erna X B 1 2 /1.76GeV) 2 ^ 'P -1 / 7 (inverting Eq. 1 of 
lBaring|[200l . where B 12 = S S /10 12 G is the scaled sur- 
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Fig. 6. — Gamma-ray luminosity above 100 MeV versus the rotational energy loss rate E. Dashed line: L 7 equal to E. Dot-dashed line: 
L 7 proportional to the square root of E. Black points: radio-loud and radio-quiet/faint pulsars. The values are extra cted from the first 
LAT pulsar catalog JAbdo et al.H2010dl) . with the addition of PSRs J0248+6021 and J2240+5832 HTheureau et al.|[20TlD . Pulsars with two 
distance estimates have two markers connected with dashed error bars. The high-_B pulsars are indicated by black squares and the values 
are reported in Table [T] 



face magnetic field and i?* is the neutron star radius. For 
PSR J1119— 6127 we observed pulsations up to e max ~ 
4GcV, which places a robust lower limit of r > 4.2 i?*, 
thereby precluding emission near the stellar surface. 

In the following subsections, we first compare the de- 
tection of PSR J1119-6127 with the two other high- 
field 7 -ray pulsars, namely PSR B1509-58 (|Abdo et al.l 
I2010ah and the rad io-quiet pulsar J0007+7307 in CTA1 
(jAbdo et all 120081) . We then will discuss possi- 
ble explanations for the non-detection in 7-rays of 
PSR J1846-0258. 

4.1. PSR Jl 119-6127 

The characteristics of PSR J1119— 6127 determined 
with Fermi are remarkably similar to those of the ma- 
jority of the ferrm-detected young pulsars (jAbdo et al.l 
l2010dl) . The delay between the radio peak and the 7-ray 
peak, the 7-ray pulse profile, and the GeV spectrum indi- 
cate an outer-magnetospherc geometry with a preference 
for the OG model as shown in the light-curve modeling 
(see M3.4L FigH]). At such high altitudes, the fields are 
far lower than B„: observe that for PSR J1119-6127, 
Bic = 5.7 x 10 3 G is somewhat high, but not as extreme 
as that of other 7-ray pulsars. This new 7-ray detec- 
tion also demonstrates that a pu lsar with higher B s than 
PSR B1509-58 can emit 7-rays (|Abdo et al.ll2010aT ) . and 
clearly shows that the emission mechanism is more com- 
plex. While PSRs B1509-58 and J1119-6127 have sim- 
ilar ages, braking indices, and surface magnetic fields, 
PSR B1509— 58 has a complex light curve with one of the 



two 7-peaks slightly preceding the radio pea k and a softer 
7-ray spectrum breaking in the MeV band (|Kuiper et al.l 
ll999tlAbdo et alJl2010aj) . Neither of these two observed 
features is fully explained by the high-altitude emission 
models. Moreover, while PSR B1509— 58 shows a non- 
thermal component up to 10 MeV and turns over in 
the MeV bands, PSR J1119-6127 presents only thermal 
emission below lOkeV but emits GeV pulsations. Com- 
paring PSR J1119-6127 with the CTA1 pulsar, it is in- 
teresting to note that apart from similarities in the GeV 
spectrum (cut-off energy at a f ew GeV) and soft X -ray 
pulsations (mainly below 2 keV. ICaraveo et"al1l2010[ ). no 
radio pul sations for CTA1 wer e detected in spite of deep 
searches (|Halpern et al.l [20041) . One likely explanation 
is that the radio beam does not sweep across our line of 
sight, and this suggests that the emission geometry (a, £) 
is different than in PSR J1119-6127. 

4.2. PSR J 1846- 0258 and other high-B pulsars 



Based on it s high E 
emission (e.g. iKuiper ~ 



= 8.1 x 10 36 erg 
Hermsedl2009T 



;s hard X-ray 
T = 1.20 ±0.01, 

20-250 keV), and its young age, PSR J1846-0258 seems 
to be a good 7-ray pulsar candidate. However, we de- 
tect no 7-ray pulsations from this object, and no source 
from the 1FGL catalog coincides with the pulsar po- 
sition. A first explanation for the non-detection of 
PSR J1846— 0258 can be simply a matter of its geometry. 
Torus modelin g in the PWN i mplies a line of sigh t angle 
C = 62° ± 5° (|Ng et all 12001 ). iLivingstone eVatl l|2009 ) 
infer a magnetic inclination a about 10° using spin-down 
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parameters applied to the model of an obliq ue rota- 
tor w ith a current-starved outer magnetosphere (|Melatosl 
1 1997T) . With these £ and a values it is unlikely that low- 
altitude radio emission will be observed. Indeed, adopt- 
ing these angles no 7-ray emission is expected from the 
OG, whereas from the TPC model non-pulsed (DC) or 
broad single-peaked pulsed emission would still be pos- 
sible (jWatters et al.l 12009?) . We stress that these pre- 
dictions are highly depe ndent on the correctness of the 
model bv lMelatos) (jl997[ ). 

Another possibility is that the sensitivity of the Fermi 
LAT is insufficient to detect PSR J1846-0258. The pul- 
sar is near the Galactic plane, 30° away from the Galactic 
center where the diffuse 7-ray background is very strong. 
According to the upper limit on pulsed flux presented in 
£ 13.51 and the distance of 7.9 ± 2.8 kpc, the derived 7-ray 
efficiency is less than 0.03 ±0.02. The Vela pulsar, which 
has a well-constrained distance around 300 pc, and shows 
a similar spin-down power (E = 7x 10 36 ergs -1 ), has a 7- 
ray efficiency of 0.01 ± 0.002, which is just at the limit of 
PSR J1846— 0258. In summary, the appreciable distance 
and high background could make the 7-ray detection of 
this pulsar difficult. 

Finally, we cannot rule out that the 7-ray spectrum 
of PSR J1846— 0258 resembles the very soft 7-ray spec- 
trum of PSR B1509-58. In X-rays, PSR J1846-0258 
mimics PSR B1509— 58 in many ways. They both 
emit pulsed non-thermal hard- X-ray emission up to 
~250 keV (~ 10 McV for PSR B1509-58) with broad 
asymmetric single-peak pulse profiles with similar hard 
spectral shapes (r ~ 1.2-1.4), and neither shows soft 
thermal emission below 2 keV (except during the 2006 
outburst of PSR J1846-0258), although such a sig- 
nal could be swamped by the non-thermal emission. 
The pulsed flux of PSR B1509-58 is - 2 x 10 -1 ° 
ergcm~ 2 s -1 (< 250 keV), which is about an order 
of magnitude higher than that of PSR J 1846-0258 
(jMarsden et al.l Il997t iCusumano et al.l 120011 ). If the 
spectrum of PSR J1846-0258 peaks in the MeV range 
as in the case of PSR B1509— 58, then the sensitivity 
of Fermi would also be insufficient to measure its pulsa- 
tions. 

The upper limits on emission from PSRs J1718— 3718 
and J1734— 3333 — which could be quiescent magne- 
tars for which we have not yet detected m a gnetic 
activity — are not constraining. I Zhang fc Jiand ()2005l ) 
studied OG ([Zhang et al.l 120041 ) predictions of energetic 
high-_B pulsars before the launch of Fermi. Assuming 
a = 55°, they predicted all pulsars in our sample to be 
detectable with the Fermi LAT. However, they assumed 
a sensitivity of 2 x 10~ 9 phcrn -2 s -1 (> 100 MeV) for 
a 2-year all-sky survey, which is not realistic along the 
Galactic plane w here the diffuse ba ckground is very 
high (see Fig. 9, lAbdo et al.l l2010dl ). Note that for 
PSR J1119— 6127 the predicted 7-ray flux is an order of 
magnitude below our observations. 

In conclusion, Fermi observations of high- 5 pulsars 
are as yet insufficient to determine whether these ob- 
jects form a homogeneous class or have both a "nor- 
mal" component and a "magnetar" component. Of the 
pulsars with fields above 10 13 G, only the soft spectrum 
of PSR B1509— 58 is anomalous, whereas the other fea- 



tures are compatible with emission from other, lower- 
field 7-ray pulsars. On the other hand, PSR J1846-0258 
has shown magnctar-like behavior, but the Fermi non- 
detection is not constraining. Note that the same high- 5 
field effects that are in place in the magnetospheres of 
these pulsars might also affect magnetars' emission. The 
non-detection of any SGR or AXP with the Fermi-LAT is 
not surprising in the context of our results. Several mag- 
netars have rotational energies and distances comparable 
to the high-5 pulsars we report her e, but no significan t 
7-ray emission has been detected yet (jAbdo et al.ll2010"cl ). 
The best hope for resolving the issue may lie with future 
X-ray polarization measurements (such as with GEM^\ 
or a successor), which can probe the spin and magnetic 
geometry of these high-i? pulsars. These measurements, 
together with deeper Fermi limits, may constrain or even 
exclude emission from the outer magnetosphere. If true, 
this would indicate that some high-_B pulsars are quiet in 
the > 100 MeV 7-ray band, and thus form a distinct class 
of objects that resemble magnetars in some of their char- 
acteristics. At present, it must also be emphasized that 
PSRs J1718-3718, J1734-3333 and PSR J1846-0258 
are distant and in the general direction of the Galactic 
center, raising the possibility that it is difficult to isolate 
their signals from the diffuse 7-ray background. 

5. SUMMARY 

We have presented the detection of PSR J1119-6127, 
which has currently the highest inferred surface mag- 
netic field of all Fermi-detected pulsars. The folded light 
curve shows a single peak which arrives with a phase 
delay of 0.43 ± 0.02 after the radio peak. No signifi- 
cant pulsations have been detected below 500 McV. The 
spectrum above 0.1 GeV can be described with a hard 
power law with an exponential cut-off (r = 1.0 ± 0.3^Q2, 

E c =0.8 ± 0.2lg;°GeV), and is similar to other nor- 
mal 7-ray pulsars, except for PSR B1509— 58. Fits of 
the pulse profile to emission models implies that the 
most likely emission geometry is that described by the 
Outer-Gap model. The best-fit parameter region (within 
radio-polarization boundary conditions) indicates a line 
of sight £ = 145° and a magnetic inclination angle 
a ~ 125°. Finally, from all these observations, the pulsar 
appears similar to a normal young 7-ray pulsar with no 
evidence for a magnctar-like behavior. 

For the other three high-S pulsars in the sample we 
have presented pulscd-cmission upper limits as a function 
of pulse width. For PSRs J1718-3718 and J1734-3333 
limits are not low enough to constrain the emission ge- 
ometry due to the low spin-down power and the distance, 
while for the excellent candidate PSR J1846-0258, the 
non-detection may suggest that the pulsar has a peculiar 
geometry, a low energy cut-off in its spectrum, or the 
Fermi LAT sensitivity is insufficient due to the distance 
and the high diffuse background. 
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TABLE 1 

Measured and derived parameters for high-Magnetic-Field Rotation-Powered Pulsars 



Parameters J1718-3718 J1734-3333 J1846-0258 J1119-6127 B1509-58 d 

Period, P (s) 3.378 1.169 0.325 0.408 0.151 

Period derivative, P (10 -12 ss" 1 ) 1.61 2.28 7.08 4.02 1.54 

Surface magnetic field, B s (10 13 G) 7.4 5.2 4.9 4.1 1.5 

Magnetic field at light cylinder, B LC (10 3 G) 0.018 0.307 13.2 5.67 42.2 

Age, t c (kyr) 33.5 8.1 0.7 1.6 1.6 

E (10 34 ergs" 1 ) 0.16 5.6 810 230 1800 

Distance, d (kpc) 4.5 ± 0.5 6.1±^ 7.9 ± 2.8 8.4 ± 0.4 5.2 ± 1.4 

Braking index, n ... 1.0±0.3 / 2.65 ± 0.01 s 2.684 ± 0.002^ 2.839 ± 0.003* 

Timing data span (months) 17 19 20 29 

Radio-7 peak offset, 5 ... ... ... 0.43 ± 0.02 0.96 ± 0.01 

7-ray peak multiplicity ... ... ... 1 2 

7-ray peak separation, A ... ... ... one peak 0.37 ±0.02 

Photon flux, a Fioo (10" 8 cm" 2 s- 1 ) < 3.8 < 3.6 < 4.0 9.3 ± 1.2 ± 2.0 

Energy flux, a Gioo (lO^ 11 erg cm- 2 s- 1 ) ... < 3.3 < 3.1 < 3.5 6.4 ±0.5 ±1.0 < 1.5 

Energy cut-off, E c (GeV) ... ... ... 0.8 ± 0.2±%% 

Spectral Index, T ... ... ... 1.0 ± 0.3±o;f 

Luminosity, b L 1 (10 34 ergs" 1 ) <8±2 < 14 ± 7 < 26 ± 19 54 ± 6 ± 8 < 4.9 ± 2.6 

Efficiency, 6 ?; < 50 ± 10 < 2.5 ±1.0 < 0.03 ± 0.02 0.23 ± 0.03 ± 0.04 < 0.003 ± 0.002 

References. - lAbdo et ail I l2010afl ; lAbdo et all d2010dT) : ^Espinoza et alj 1 1201(H) ; ^Livingstone et al.l 1 120061) ; ' iWeltevrede 
ILivingstone et al.l ( 1200511 

Note. — The upper limits on pulsed flux for PSRs J1718-3718, J1734-3333, and J1846-0258 assume a Gaussian peak with a FWHM of 

profi le, and an exponent ially cut-off power-law spectrum (sec Eq. 0) with T — 1.5, E c — 3 GeV, and /3 — 1. The upper limit for PSR. B1509 
from lAbdo et al.l I l2010al) using the LAT 0.1-0.3 GeV upper limit in Fig. 3. 

a E >0.1 GeV 

fn is assumed to be 1 which can result in an efficiency > 1. 
c The errors on the upper limits are dominated by the distance uncertainties. 



J0007+7303 e 



0.316 
0.36 
1.1 
3.21 
13.9 
45 

1.4 ± 0.3 



0.23 ±0.01 

30.7 ± 1.3 
38.2 ± 1.3 
4.6 ± 0.4 

1.38 ± 0.05 
8.9 ± 3.8 

0.20 ± 0.08 
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